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The application of spatially- uniform light on conventional insulators can induce Floquet spectra 
with characteristics akin to those of topological insulators. We demonstrate that spatial modulation 
of light allows for remarkable control of the properties in these systems. We provide configurations to 
generate one-dimensional bulk modes, photoinduced currents, as well as fractionalized excitations. 
We show a close analogy to p-wave superconductors and use this analogy to explain our results. 



The Quantum Hall Effect [lj lead to the discovery of 
a close connection between topology and certain physi- 
cal properties of condensed matter systems [2H4]. Our 
understanding of the role of topology has greatly ex- 
panded following the recent discovery of new classes of 
topological phases and of new materials displaying topo- 
logical properties [5HT2]. Topological phases are char- 
acterized by integer-valued numbers that are invariant 
to small changes of their Hamiltonian. This makes in- 
triguing effects, such as quantized Hall conductivity and 
non-abelian excitations, robust properties of these sys- 
tems HJ HMD- 

Recently, it has been shown that topological properties 
can be induced in conventional insulators by the applica- 
tion of time-periodic perturbations fT6ll22] . Proposals for 
these so-called "Floquet Topological Insulators" (FTIs) 
include a wide range of physical solid state and atomic 
realizations, driven both at resonance and off-resonance. 
These systems display metallic conduction enabled by 
quasi-stationary states at the edges fT6| \T7\ [20] . Dirac 
cones in three dimensional systems [21], and even Flo- 
quet Majorana fermions [22 J. 

In this Letter, we demonstrate dramatic effects that 
arise in FTIs when light is modulated in space. Non- 
uniform light can give rise to controlled one-dimensional 
modes in the bulk, to fractionalized excitations, and to 
photoinduced electric currents. We establish these re- 
sults both numerically and analytically. We show that 
the Floquet spectrum resembles that of a p-wave super- 
conductor with a spatially- modulated order parameter. 
This analogy provides a simple description of the mech- 
anism behind our results. We propose setups by which 
the properties of light-induced topological phases can be 
controlled. For example, by modifying the angle of in- 
cident light on a system one can set the density of one- 
dimensional modes in its bulk. 

We begin by building a description of FTIs in a generic 
zincblende lattice model. The unperturbed system is 
given by the Bloch Hamiltonian 



H k 
HI 



(1) 



This can describe, for example, HgTe quantum wells, 



in which case Hi? 



is a 2 x 2 Hamiltonian 



acting on the subspace spanned by the J z = (^, §) 
[Jz = (— §)] states, respectively. Thus, the two 
blocks in Eq. (HI) are related to each other by time re- 
versal symmetry. Most generally, one can write 



Hk = dk • <r + £khx2- 



(2) 



We consider time-dependent perturbations that do not 
connect the two Hamiltonian blocks, and perform the 
analysis on a single block. For example, we will study 
the 2x2 Hamiltonian 

Hu n (t) = d k • <7 + Ekhx2 + Vk ■ <rcos (ujt + a) . (3) 

where Vk can depict the effect of shining linearly polar- 
ized light on the sample [20J. 

The solutions of the Schrodinger equation for a 
time-dependent system evolve according to i/)(t') = 
U(t f ,t)ip(t), where U is the time evolution operator 



U(t',t) 



exp 



H (t") dt" 



(4) 



For a time-periodic system, Floquet 's theorem states that 
these solutions can be written as ip (t) = e l£at cp a (t) , 
where (p a (t) = cp a (t + r) and r = ^ [23J. The e a are 
called quasi-energies, and are only denned modulo u; the 
ip a (t) satisfy the eigenvalue problem HF(p a (t) = s a (p a {i), 
where Hp is the Floquet Hamiltonian, obtained from the 
time evolution operator over a full cycle 



-iHft 



U(t,0). 



(5) 
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Figure 1: fife, defined in Eq. (pj, for the lattice model, Eqs. 
Q and ([9]), with a = 0. (a) x and y components of nu 
. (b) n\. Note that Hk is in a hedgehog configuration, as 
it wraps the unit sphere exactly once. This corresponds to 
Cf = 1. The dashed lines depict the range [— ko,ko] over 
which C' kx ^ 0. 



Figure 2: (a) Floquet spectrum of Eq. ^ for a domain wall 
configuration. The plotted range corresponds approximately 
to [— ko,ko].(b) Amplitude of the localized bulk (red) and 
edge (blue) states . Inset: Floquet spectrum for constant a. 
Results are for V& = A = -B = 0.2M, uj 2.7 M. 

For a 2 x 2 block, Hp can be written most generally as 
H F = n k ■v + Skhx2- (6) 

This has the structure of a gapped system provided that 
n k does not vanish on the Brillouin zone. We then intro- 
duce a topological invariant [7J [5] 

C f = ^~ J J ( d k x n k x d ky n k ) • n k (7) 

where h k = n k /\n k \. Cp can be nontrivial even when 
the unperturbed system is topologically trivial [20J . As a 
consequence, the time-dependent perturbation can give 
rise to topologically-protected edge states in Hp. 

The Floquet spectrum is independent of the value of 
a G [0,27r]. To see this, note that a shift in a results 
in U a (r,0) = W^U a=0 (r,0)W : where W = ^ (g,0). 
Thus, changing a is equivalent to a similarity transfor- 
mation of H F , H F (a) = W^H F (0)W. The spectrum of 
Hp is therefore independent of a, even if its eigenstates 
are not. In particular, varying a cannot induce a gap 
closure, nor change Cp. Naively, this implies that at the 
interface between two portions of the sample with differ- 
ent values of a, no edge modes are expected in the Flo- 
quet spectrum. However, we'll show below configurations 
in which localized modes appear at such an interface. At 
first glance this is remarkable, since the interface connects 
two systems with identical topological classification. 

Lattice model- We first demonstrate the existence of 
quasi-stationary interface modes in a lattice model. Con- 
sider the time-independent Hamiltonian, Eq. (J2|, with 

d k = (A sin k x , As'mk y , M + 2B (cosk x + cosk y — 2)) . 

(8) 

We choose e k = 0. Then Eq. Q has particle-hole sym- 
metry (PHS) in which the valence and conduction bands 
are interchanged. We add a time-dependent perturbation 

V = V a z cos (ut + a) . (9) 

We evaluate the time evolution operator numerically, 
Eq. Q, by discretizing the time interval between t = 
and t' — r = — . We then obtain the Floquet spectrum 



Figure 3: Amplitude of the vortex bound state. Inset: quasi- 
energy of the state as a function of the system size L. Results 
are for V 1.5M, A 0.2M, B -0.1 M, uj 2.7 M 

by diagonalizing U, see Eq. When a is a constant, 
the system has translational invariance and n k can be 
calculated for each k separately. Fig. [l] shows that h k 
wraps exactly one time around the unit sphere over the 
first Brillouin zone. The system is therefore topologically 
non-trivial with Cp = 1, and we expect one localized chi- 
ral mode at each edge of the system. To demonstrate this, 
we choose cylindrical geometry with open boundaries in 
the y direction and periodic boundaries in the x direction. 
Figure [2] shows that indeed, one zero quasi-energy mode 
exists at each boundary of the system. This corresponds 
to the edge state found in Ref. [20] , 

We now allow a to be position-dependent. As a first 
example, we consider a domain wall configuration, across 
which the external perturbation changes sign, a (y) = 
7T0(y). This captures the phase shift across the nodes 
of a standing wave created by two interfering light rays 
incident on the sample. By adjusting the incidence angle 
of the rays, one can control the periodicity of the stand- 
ing wave such that these nodes are well-separated. We 
choose cylindrical geometry. Since the system remains 
translationally invariant along x, we work in the hybrid 
coordinate basis (k x ,y) and diagonalize U for each k x 
value. Figure [2] shows the resulting quasi-energy spec- 
trum. Note that, in addition to the edge modes, the spec- 
trum now includes two counter-propagating zero quasi- 
energy modes localized near the domain wall at y = 0. 

As a second example, we consider a vortex configura- 
tion, in which the phase a winds by 2tt about a point, 
a(f) = arctan(y/x). A lattice of such vortices can be 
created by interfering three lasers, and an isolated vor- 
tex can be created using a phase mask. We set open 
boundaries and diagonalize U in real space to obtain the 
Floquet spectrum. In addition to the edge modes, we 
find a zero quasi-energy state bound to the vortex. Fig- 
ure [3] shows the wave function of the bound state. This 
state hybridizes with an edge state and opens a small gap. 
However, this is a finite-size effect, and the gap energy 
decays exponentially with system size. 

Analogy to p x + ip y superconductors. - In order to ex- 
plain these results, we establish an analogy between our 
system and a p x + ip y superconductor (pSC). We then 
relate the domain wall and vortex core states to the well- 



3 



known zero energy states of pSCs. 

We derive an approximate analytic expression for Hp. 
The derivation is first carried out for constant a. We 
start by omitting the component of Vk parallel to d^, 
since it only affects the dynamics weakly when averaged 
over a full cycle and is known not to influence the topo- 
logical properties of the system [20]. The remaining per- 
pendicular component is Vk ± = Vk — {V k -dk)d k . We 



We apply these results to our system. The unper- 
turbed Hamiltonian for small k is Hk = Ak • a ± Ma z . 
The Floquet Hamiltonian in the RWA is then 



define Vk 



and introduce Vk = Vk t ±/Vk t ± and 



Wk — dkXVk, such that dk, Vk and Wk form a right-handed 
triad. The perturbation can be decomposed into terms 
that rotate and counter-rotate about <4> 

Hk ~ dk • <r + \Vk,±_ (vk cos out ± Wk sin out) • a 
+ |T4 ? _l (vk cos out — Wk sin out) • a. 



We go to a rotating frame through the time-dependent 
unitary transformation R(t) = exp (— idk • • The 
resulting states \ip(t)) r = satisfy id t \^(t)) r = 

# 7 |^(t)) r , where H' = (il - ild^j R. We then find 

H k = (d k - %d k ± |Vfe,_L (fife cos a + sin a)) • a 
+ \Vk,±_ (vk cos (2out + a) — Wk sin (2out ± a)) • a. 

The Rotating Wave Approximation (RWA) consists of 
neglecting the 2ou term in H' k . This is valid near a res- 
onance, \ou — A£?| <C ou ± A£?, where A£? is the energy 
difference between states in the lower and upper bands. 
This procedure yields a time-independent operator 



ttRWA 
H F,a 



(( 1 ~ 2^) dk + \Vk t ±. (i)kCosa + w k sin a)^j • a, 

(10) 

which is the Floquet Hamiltonian in the RWA. 

The analogy to a pSC can be seen explicitly by per- 
forming a /^-dependent unitary transformation that ro- 
tates (dk,Vk,Wk) — » (z/k,z x k), where fc = (k x ,k y ). 
This leads to a Hamiltonian of the form 



(fee 



A fc e" ia - ik y ) 
■ ik y ) -Ck 



(11) 



where (k = dk — ou/2 and A& = Vk t ±/2k are real. Equa- 



tion (11) resembles the Hamiltonian of a pSC with com- 
plex order parameter Ake~ za . The analogy to the pSC 
can be seen graphically in Fig. [TJ where panel (a) depicts 
the normal dispersion and panel (b) the superconducting 
order parameter, which is seen to have p x + ip y symme- 
try. Note that, unlike an actual superconductor in which 
the Nambu basis describes particle and hole states, here 



Eq. (11) acts on two particle-like states corresponding 



to valence and conduction bands of the Floquet problem. 



Hence, the spectrum of Eq. ( 11 ) matches the correspond- 



ing pSC, but the nature of the wave functions in the two 
cases is related by a particle-hole transformation. 



ttRWA 
H F,a 



(jj (z + jtfk^j + A (kcosa + z x /csina 



where rj = (M — ^) and Aq = 4^-. When a varies 
in space, this is generalized to a Bogoliubov-de Gennes 
(BdG) equation (here ip = (u,v) T ): 



(e-rj)u= (±r) - A (r)) * (-id x + d y ) v 

A* (f)) u 



(e + rj) v = (-id x 



(12) 



where * denotes the symmetric product a * b = 
7}(ab + ba). In theories of superconductivity the gap 
function A (r) is calculated self-consistently. Here, A (r) 
is directly determined by the external perturbation. 

When a (r) describes a vortex, the BdG equation has 
a zero energy solution bound to the vortex core. This 
state is analogous to the well-known pSC vortex core 
modes, and just as in the case of a pSC, it is topo- 
logically protected provided particle-hole symmetry is 
present [24, 25J. This gives a fractional charge of ±1/2 at 
the vortex, depending on whether the state is occupied 
or not [26]. In the full 4x4 system, the vortex is felt in 
both subspaces, yielding two zero energy states which can 
be independently occupied. This leads to fractionalized 
excitations (±1/2, ±1/2), where the indices denote the 
charge in each subspace. The scenario is closely anal- 
ogous to fractionalization in polyacetylene, where frac- 
tional charge excitations in the spin-up and spin-down 
channels combine into spin 1/2 excitations with integer 
charge [27]. 

When a (f) describes a domain wall, the BdG equation 
has two zero quasi-energy solutions, localized around y = 
with a localization length of £ = . This configuration 
is analogous to an interface between two pSCs with a 
relative phase difference of ir - a u tt junction" - which 
is known to have zero energy bound Andreev states [25j 
28J. To understand these modes, let's first consider a 
system with no disorder and no interactions. Then k x 
is a good quantum number and the system reduces to a 
ID pSC for each k x . We can then define a k x - dependent 
topological invariant, C k ^ as the winding number of 
in the (n k: n k ) plane. As can be seen from Fig. [lj when 
a = the ID pSC is topological with C k = 1 for k x G 
[— fco? ko] and it is trivial otherwise. On the other side of 
the domain wall, a = 7r, and the sign of C kx is reversed. 
For k x G [— fco? ko] the change in sign in C kx implies a pair 
of localized states exists at the interface, which disperse 
in opposite directions as a function of k x . In particular, at 
k x = these modes cross with zero quasi-energy. When 
PHS is present, one of these states is even under the 
PH transformation, whereas the other is odd, preventing 
them from mixing and opening up a gap. This protection 
is robust even in the presence of disorder and interactions, 
provided these do not break PHS. 
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Note that the domain wall and vortex configurations 
are closely related. To see this, we write a (r) = 
arctan(y//te), which describes a vortex for /? = 1 and 
a domain wall for /3 — >• + . Thus, the domain wall is 
a continuous deformation of a vortex and the ir jump 
arises from the vortex winding. The interface modes are 
therefore smoothly connected to the vortex core states. 

Photoinduced current - Recall that imposing a phase 
twist on a superconductor Ae m ( r ) results in a Joseph- 
son current js — p s ^ a [31] • Motivated by the analogy 
to pSC, we consider the effect of a slowly- varying phase 
twist a(y) = + yd y a, where d y a is a small constant. 
This can be achieved, for example, by shining a coherent 
light ray incident at an angle to the surface of the sample. 
Indeed, we find that the system experiences a DC cur- 
rent along y, and compute an analogue of the superfluid 
stiffness p s in terms of the vectors dk and n^. 

By Noether's theorem, the current density operator is 



where V& = (dk x ,dk y )- Our goal is to compute the ex- 
pectation value of Eq. ( [l3| ) with respect to the state at 
half-filling, in which all eigenstates in the Floquet valence 
band are fully occupied, 



(14) 



where the sum is over all the negative quasi-energy states. 
Here, we insert R(t) — exp (—idk • cr^j since the cur- 
rent is computed in the lab frame. 

In order to derive an expression for the current, we 
compute Hp to linear order in d y a in the RWA. We ob- 



tain Hp = H 



by Eq. (10) and 



RWA 

F,a 



(d y a)H? WA , where H^a ig given 



RWA 



y 



Vk,± (cosaoWk — sinao^/e) • ~c? + h.c. 



We then write \ipk)r = Wk)r + 9 y a\^l) r , where |^) r 
are eigenstates of Hp^ Q A , and \ip\) r is obtained from 
I s * order perturbation theory. The resulting current, to 
O (d y a), vanishes along x, while along y it is: 



i = ii <gfr& (d k x vj k ) ■ v k r (r k \ iV m 



where d = dip) — (dip\ ip). In this expres- 
sion, we have summed over the contributions coming 
from both sub-Hamiltonians, and H*L k , appearing in 
Eq. In terms of dk and n^, this can be rewritten as 



(jy) ff d 2 k Vp 
d y a JJ (2tt) 2 2 



l^(^V^) 2 ( 15 ) 



where di- = 



Integrating Eq. (15) over the Brillouin zone gives a 
non-zero DC current. As a check of our analysis, we 
compared these results with numerical simulations on the 
lattice model and found good agreement. This result is 
reminiscent of the photogalvanic effect proposed to exist 
at the surface of 3D topological insulators radiated with 
circularly polarized light [32J. 

Discussion.- The analogy to pSC relies on particle- 
hole symmetry. Thus, it is natural to ask which of our 
results rely on this symmetry. For instance, Eq. (15) 



for the photoinduced current is valid even when PHS is 
broken. In contrast, PHS plays an important role in pre- 
venting the interface modes from opening a gap, as dis- 
cussed above. For real systems PHS is only approximate, 
and we must consider the effects of weak PHS breaking. 
We find that this induces a small mixing of the interface 
modes. For example, using the same parameters as in 
Fig. [2j but adding a PHS breaking term to Eq. ^ of the 



form £k 



-0.2 (cos k x + cos k y — 2), we find that a small 



and n k Ind 



(13) gap w 10 -5 opens up. Thus, experiments carried out at 
temperatures above this gap will not be sensitive to PHS 
breaking. The vortex core state shows higher degree of 
robustness to breaking of PH symmetry. Although PH 
symmetry breaking shifts the bound state quasi-energy 
away from zero, it remains a mid-gap state provided the 
symmetry is weakly broken. In particular, the fractional 
nature of the excitation is unaffected. 

We now briefly discuss the physical manifestations 
of our results in realistic systems. Earlier work has 
shown that resonantly-driven systems can reach a steady 
state, with occupation given by the Fermi-Dirac distribu- 
tion based on the Floquet spectrum, and with an effec- 
tive temperature that is determined by interaction and 
phonon relaxation rates [331 El] • Then, a standing wave 
pattern is expected to create one dimensional channels 
along the domain walls and thus induce anisotropic con- 
ductivity. It may be possible to generalize our results to 
systems in which the light frequency exceeds the band- 
width (that is, systems that are driven off-resonance). 
The advantage of doing this is that physical properties 
of systems driven off-resonance are much easier to under- 
stand [T^ \T7\ 155], Modulated FTIs driven off-resonance, 
and in three dimensions, will be discussed in future work. 

In conclusion, we provided a set of schemes by which 
the properties of Floquet topological insulators can be 
manipulated using modulated light. We proposed ex- 
plicit setups by which bulk ID channels, fractionalized 
excitations, and light induced currents can be generated. 
Our analysis demonstrates the great potential for tun- 
ability and control of light-induced topological phases. 
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